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Monte-Carlo study of the Ap parameter* 
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We present results concerning a lattice study of the electroweak p-parameter. We have used an SU(2) x U(l) 
symmetric chiral Yukawa model built with Zaragoza fermions. The decoupling of the species doublers in this 
model is verified numerically. We find that the numerical data for Ap are well described by one-loop perturbation 
theory in the same finite volume and with the same finite cut-off. However, a finite cut-off can cause substantial 
deviations of Ap from the standard value, even in infinite volume. 



1. INTRODUCTION 

An important quantity in the Standard Model 
of electroweak interactions is the p-parameter, de- 
fined as the relative strength between the neu- 
tral and charged current interactions at zero mo- 
mentum transfer. At tree level at low momen- 
tum, the electroweak interaction is described by 
the standard current-current interaction and it is 
well known that pt iee = M^/M^ cos 2 0\\r = 1- 
When quantum effects are taken into account, 
p = (1 — Ap)" 1 deviates from 1. We will be 
concerned here only with the dominant "univer- 
sal" (i.e. process-independent) part of Ap coming 
from the difference in vacuum polarization of the 
W and Z propagators. 

In the Standard Model, the dominant contri- 
bution comes from the t-b mass splitting. In the 
limit m 2 <C m 2 it is given by 
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From a phenomenological point of view, this 
quantity is an interesting tool to investigate the 
existence of unknown physics through loop ef- 
fects. This is exemplified by the bounds put on 
the top quark mass before its discovery. As it 
appears in eq. (0), there is a strong dependence 
on the mass of a heavy fermion running in the 
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loop. In particular such a fermion does not de- 
couple from the low-energy physics. It would be 
desirable to control the persistence or not of this 
non-decoupling phenomenon non-perturbatively, 
for large values of the mass, and to study its re- 
lation with the issue of triviality. This problem 
has been addressed in a large- Np approximation 
before §. 

From the lattice point of view, these character- 
istics have some relevance to the doubling prob- 
lem. We know for example that heavy doublers 
treated a la Wilson- Yukawa do not decouple in 
the continuum limit. Through the relation ([!]), 
Ap can be used as a counter of fermions. It thus 
provides us with a tool to test models of chiral 
lattice fermions. 



2. MEASURING Ap ON THE LATTICE 

A numerical prediction for Ap from a direct 
study of the gauge boson propagators would re- 
quire a tractable lattice formulation of a four- 
dimensional chiral gauge theory. This task is be- 
yond our present capabilities. However, the gauge 
interactions are weak compared to the Yukawa 
and scalar ones and it is a reasonable approxi- 
mation to neglect all the gauge couplings. There 
are Ward identities relating the vacuum polar- 
ization tensors of the gauge bosons to those of 
the charged and neutral goldstone particles in the 
fermion-scalar sector, and p can be expressed as 
the ratio of their wave function renormalization 
constants: p = Zq/Z + Bpf. 
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In the gaugeless approximation, therefore, only 
the scalar-fermion sector of the Standard Model 
needs to be latticized to obtain p. For this pur- 
pose, we have used an SU(2) x U(l) symmetric 
chiral Yukawa fermion model built with Zaragoza 
fermions. The lattice action is: 
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$ is a fixed modulus Higgs field ($ G SU(2)), V> is 
a fermionic field for two fermion doublets, Pr is 
the right projector 5(1+75) and Y — ( v ° ) is the 
matrix of Yukawa couplings, where y = yt op and 
2/bottom is set to zero. The fermion field ip^(x) is 
the average of the ip fields over the vertices of the 
hypercube at a; [[| ; in momentum space it can be 
written as i/j^Hp) = ■P'(p)V'(p) with "form factor" 
F(p) = JT A cos(apA/2). This model is expected to 
describe two coupled fermion doublets and thirty 
massless and uncoupled doubler doublets . We 
have simulated the action (||) with the Hybrid 
Monte-Carlo algorithm (the reason why we have 
doubled the fermion content) and present results 
for 6 3 x 12 and 8 3 x 16 lattices. Details will be 
given in Q. 

3. DECOUPLING OF THE DOUBLERS 

With Zarogoza fermions, the decoupling of the 
doublers is known to hold in perturbation theory 
PJ|. Numerically, the structure of the phase di- 
agram was previously found in good agreement 
with a mean-field prediction, giving some indica- 
tion that the decoupling works (j| . We have done 
some further tests reported here: 
• From the three point function Gz{k,p, q) = 
Ei.s.z^ot 1 )^!/)^ 2 )) exp[-i(k.x + p.y - q.z)], 
where IIo(x) = Tr [$(x)cra] is an interpolating 
field for the neutral goldstone, we have extracted 
the renormalized three-point couplings ya p t for 
the usual fermion (i.e. with momentum near 
p = 0), and ydb for a doubler (i.e. with momen- 
tum near p^ = = na^ 1 ) respectively. We find 
as expected that j/db ~ 2/3pt *F 2 , i.e. the coupling 
of the doublers is suppressed by the form factor 
F 2 . This suppression factor is around 10 -2 for 
the minimum momentum allowed for a fermion 



on an 8 3 x 16 lattice and will decrease if the vol- 
ume is increased. 

• We have studied the one-loop renormalization 
group equations for the evolution of the scalar 
and Yukawa coupling constants for the contin- 
uum model with 2 and 32 fermions doublets and 
compare the solutions with our Monte-Carlo data 
in fig. [|. The data seem to follow the one-loop 
prediction and support the expectation that we 
have only two dynamical doublets, without cou- 
pled doublers. 
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Figure 1 . Higgs mass as a function of the fermion 
mass in units of v r . MC data for 6 3 x 12 and 
8 3 x 16 lattices are shown, and the solid line is 
the perturbative result with 2 doublets, for a scale 
ratio log(i; r /A) = —4/3. The analogous result for 
32 doublets is given by the dotted curves. 



• We have compared results obtained for Ap 
with Zaragoza and naive fermions for two points 
in the phase diagrams with comparable fermion 
mass in lattice units. We defined the renormal- 
ized coupling constant, y r , as the ratio of the 
renormalized fermion mass divided by the renor- 
malized vacuum expectation value of the scalar 
field, v r . On a 6 3 x 12 lattice, the compari- 
son between Monte-Carlo data and perturbation 
theory gives [Ap MC /y 2 ] / [A Ppcrt /y 2 ] = 0.98(13) 
for k — 0.24, y=0.6, am/=0.311 in the naive 
case and [Ap MC /y 2 ] / [Ap pcrt /y 2 ] = 0.92(18) for 
K = 0.31, y—1.0, etm/=0.345 in the Zaragoza 
case. This means that for both fermion mod- 
els, the numerical results are well compatible with 
their respective perturbative predictions. On the 
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other hand, perturbation theory predicts that 
A/Onaive/Apzaragoza equals 17.13 and 17.26 for 
arrif equal to 0.311 and 0.345 respectively (this 
ratio becoming 16 in the infinite volume limit 
when arrif <C 1). This is a clear indication that 
the Zaragoza doublers are decoupled and do not 
contribute to Ap. 

4. RESULTS FOR Ap 

To detect possible non-perturbative effects, we 
have compared the results obtained from the sim- 
ulations with the one-loop perturbative predic- 
tions at finite volume V and finite cut-off 
a -1 . Namely, for each lattice (re, y, V) we have 
extracted from the numerical data the value of 
the fermion mass in lattice units, arrif, and used 
it in a finite- volume fermionic one-loop perturba- 
tive computation of Ap. Fig. || shows the ratio of 
the Monte-Carlo value ApMc/lJr divided by the 
corresponding perturbative prediction. The nu- 
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5. CONCLUSIONS 

The dominant part of Ap can be computed 
within the scalar-fermion sector of the Standard 
Model. This can be done on the lattice with 
present techniques. We have found different kinds 
of numerical evidence that the decoupling mech- 
anism works as expected with the Zaragoza pro- 
posal for lattice fermions. Ap is well described by 
the perturbative prediction. Nevertheless, even if 
Ap is always given by perturbation theory, cut- 
off effects can be important [Q. When the mass 
of the fermion approaches the cut-off, Ap devi- 
ates from the standard result given in eq. (Q), 
valid when m//A <C 1. In fig. [| we give the ratio 
Ap peT t(airif) I A/? per t(0) for the lattice model we 
have used here. It shows that cut-off effects in 
Ap can be quite substantial. 
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Figure 3. Finite-cut-off effects in Ap. We plot 
the perturbative one-loop estimate for the ratio 
Ap(arrif)/ Ap(0) obtained with the action (2) at 
infinite volume. 



Figure 2. Ratio of the Monte-Carlo Apjyl di- 
vided by the finite-volume, finite-cut-off pertur- 
bative value. Data from 6 3 x 12 and 8 3 x 16 lat- 
tices are shown. 



mcrical results seem to be well represented by the 
finite-volume, finite-cut-off one-loop perturbative 
calculation for both volumes and all the values of 
y so far explored. 



A more detailed account of this work, including 
a more complete list of references and acknowl- 
edgements, will appear soon |J. 
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